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We propose a model of suppression of 7r°'s based on two different effects: at low pr we take 
into account the shadowing corrections, while at large pr the suppression is produced due to 
the interaction of the large pr pion with the dense medium created in the collision. The main 
features of the data on AuAu and dAu collisions are reproduced both at mid and at forward 
rapidities. 



1 Introduction 

The experimental data on centvalAuAu collisions taken at RHIC show that the nuclear modifica- 
tion factor, RAA(b,y,pT) = dyd'2 pT W / \_ n (b) dyd^p T > * s smaner than 1 for all px- That means 
that the yield of particles produced in AA collisions at mid-rapidities and large pt increases 
with centrality much slower than the number of binary collisions n(b). This phenomenon is 
particularly interesting since it is not observed in dAu collisions at RHIC at mid-rapidities. 

In order to explain these data, one should take into account two kinds of effects. On one side, 
at low px this suppression is explained by the shadowing corrections^. At very high energies 
the shadowing effects -nonlinear- lead, as s — > oo, to a saturation of the distributions of partons 
in the colliding nuclei. These corrections are important for the description of inclusive spectra, 
for particles with px ~< Pt >• But with increasing px the shadowing corrections decrease and 
the scaling with nib) is predicted in perturbative QCD (Raa — ► 1)- 

So in order to explain the data at large px, it is neccessary to consider that for particles 
with large momentum transfer, there are in general final state interactions^ -jet quenching or 
jet absorption-. These interactions lead to an energy loss of the large px parton (particle) in 
the dense medium produced in the collision. Hadrons loose a finite fraction of their longitudinal 
momentum due to multiple scattering. From this point of view, it is natural to expect that a 
particle or a parton scattered at some non zero angle will also loose a fraction of its transverse 
momentum due to final state interactions with a medium of partonic or prehadronic nature. 
The scattered particle does not disappear as a result of the interaction but its pt is shifted to 
smaller values. Due to the steepness of the px distribution, the effect of suppression at large 
px may be quite large. Moreover, in this case there is also a gain of the yield at small px due 
to particles produced at larger px - which have experienced a px shift due to the interaction 
with the medium. This suppression vanishes at low px- when the px of the particle is close 



to < pt >, its pt can either increase or decrease as a result of the interaction, i.e. in average 
the pt shift tends to zero. Our aim here is to describe the suppression of the yield of pions in 
a framework based on final state interactions and taking into account shadowing and Cronin 
effect. 



2 The model 



Our approach contains dynamical, non linear shadowing. It is determined in terms of diffractive 
cross sections and it applies to soft and hard processes. The reduction of multplicity from 
shadowing corrections can be expressed by 

Jdhf A (s)f B (b-s) T A {b) 
RAB(b) = TXjM ' Uh) =l + AF{s)T A {b) (1) 

where the function F(s) represents the integral of the triple Pomeron cross section over the 
single Pomeron one: 
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F(s) = 4vr / dy ™ ° | t=0 = C[exp(y max ) - exp(y min )]. (2) 

From this formula one can observe that for a particle produced at y = 0, with y ma x = 5 m (s/my) 
and y m in = ln(R A m,N /V~3), the effect of shadowing decreases at large pr due to the increase of 
m T in y ma x- 

In other words, due to coherence conditions, shadowing effects for partons take place at 
very small x, x <C x cr = I/ttinRa where mjv is the nucleon mass and R A is the radius of the 
nucleus. Partons which produce a state with transverse mass tut and a given value of Feynman 
xf, have x = x± = \{\J x 2 F + 4m T /s ± xf)- Our shadowing applies to soft and hard processes. 
Nevertheless, for large pt these effects are important only at very high energies, when x ~ ^t- 
satisfies the above condition. At fixed initial energy (s) the condition for existence of shadowing 
will not be satisfied at large transverse momenta: In the central rapidity region (y* = 0) at RHIC 
and for pr of jets (particles) above 5(2) GeV/c the condition for shadowing is not satisfied and 
these effects are absent. 

So for the large pt partons it is neccessary to include the final state interactions. The 
interaction of a large px particle with the medium is described by the gain and loss differential 
equations which govern final state interactions: 

dp H (x,p T ) ~ r / \ / i_ jf \i a\ 
= -a p s [ph(x,Pt) ~ Ph{x,Pt + opT)\ (3) 

where ps and pu correspond to the space-time density of the medium and of large pr vr°'s, a is 
the interaction cross-section averaged over momenta and d 4 x represent the cylindrical space-time 
variables: longitudinal proper time r, space-time rapidity y, and transverse coordinate s. 

Assuming a decrease of densities with proper time 1/r -isoentropic longitudinal expansion, 
transverse expansion is neglected- the above equation tranforms into: 

r dN «°( b £y>Pr) = -aN(b,s,y) [N n oQ>, s,y,pr) ~ N w »(b, s,y,p T + 8pr)\ (4) 

where N(b, s, y) = dN/dy d 2 s(y, b) is the density of the medium per unit rapidity and per unit 
of transverse area at fixed impact parameter, integrated over pr, and N n o(b, s, y,pr) is the same 
quantity for 7r°'s at fixed px- 

If we integrate from initial time ro to freeze-out time r/ and taking into account the inverse 
proportionality between proper time and densities, tj/tq = N(b,s,y)/N pp (y), where N(b,s,y) 



is the density produced in the primary collisions in DPM and N pp (y) is the density per unit 
rapidity for minimum bias pp collisions at yfs = 200 GeV= 2.24 fm~ 2 , we obtain the suppression 
factor S n o(b,y,px) of the yield of 7r°'s at given px and at each impact parameter, due to its 
interaction with the dense medium. We get: 

c n J d 2 s a AB {b) n(b,s) S w o(b,s,y,p T ) 

J d z s a AB (b) n(b,s) 

where the survival probability is given by: 

N n o(b,s,y,p T + 5p T ) 



S w o{b,s,y,p T ) = exp <^ -a 
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N w o(b,s,y,p T ) 



Here a AB (b) = {1 - exp[-a pp AB T AB (b)]}, T AB (b) = J d 2 sT A {s)T B (b - s), T A {b) are the profile 
functions and n(b, s) = AB a pp T A (s) T B (b — s)/a AB (b). 

One can estimate which amount of the effect takes place in the partonic phase and which 
one happens in the hadronic phase. We can divide our suppression factor in two parts: 
Partonic, from initial density N(b,s,y) = ~ I2j» to ^10L ~ or equivalently from 

To = 1 fm to t p = 3.36 fm, and hadronic, from partonic density ^jgr ~ ^p- to N pp (y) = 

d ^^ y = 2.24 fm -2 , or equivalently from t p = 3.36 fm to rj = 5 — 7 fm. We find that 75% of the 

effect takes place in the partonic phase while only 25% of the effect takes place in the hadronic 
phase. 



3 Numerical results 



In order to perform numerical calculations, we need a -our free parameter, a ~ 1 mb- and the 
PT distribution of the 7r°'s. We have proceeded as follows: 

In pp collisions at = 200 GeV, the shape of the px distribution of 7r°'s can be described 
as (1 + PT/po)~ n with n = 9.99 and po = 1.219 GeV/c. The corresponding average pr is 

< Pt >= 2po/( n — 3) = 0.349 GeV/c. In AuAu collisions we assume that the px distribution 
of 7r°'s at each b is given by the same shape, (1 +Pt/po)~ u , with n = 9.99 and changing the 
scale po into po(b) =< pt >b ( n — 3)/2, where < px >b is the average value of px measured 
experimentally at the corresponding centrality. For central AuAu collisions {n part = 350 ,6 = 2), 

< Pt >exp is equal to 0.453 GeV/c resulting in a po of 1.583 GeV/c. 

We take n as fixed, since at large px the shadowing vanishes and the ratio R AA is independent 
of px- Also, the experimental value of n in dAu is the same as in pp within errors. 

We can thus compute the ratio R AA in the absence of final state interaction (colunm I, Table 
1). In this case R AA increases with px, there is shadowing are low px and we have a Cronin 
effect that does not decrease at large px- 

To these values we apply the correction due to the suppression factor S n o: 

First, we neglect the second term -the gain- in eqs. (3-4) (colunm II, Table 1). We find a 
general suppression in R AA , independent of px- Raa increases slightly with px and agrees with 
experimental data for px > 5 GeV/c. At lower px the result is significantly lower than the data. 
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Table 1. Values of R\ uAu {p T ) for the 10 % 
most central collisions AuAu collisions at mid- 
rapidities (|y*| < 0.35). Column I is the result ob- 
tained with no final state interaction. The results 
in the other columns include final state interac- 
tion with several ansatzs for the px shift induced 
by this interaction. 



Then we introduce the second term -the gain- in eqs. (3-4) in two ways: 
Assuming that the px shift of the 7r°, due to its interaction with the medium, is constant 
(two cases : 5px = 0.5 GeV/c and 5pr = 1.5 GeV/c) (columns III and IV), we obtain a slight 
increase of Raa at large px, rather insensitive to the value of the shift, consistent with data. 
The problem at small px remains. 

If we assume that 5pT oc {px~ < Pt >b) (columns V and VI), in such a way that the factor 
S n o is 1 at pt =< Pt > as it should be -no suppression at low px~, we obtain a slight decrease 
of Raa at large px, consistent with data and an increase of Raa at low px -due to the shift of 
large pt particles- resulting in an agreement with data. 
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Figure 1. Left: Values of Ra u Au(pt) for the 10 % most central collisions (lower line) and for peripheral 
(80-92 %) collisions (upper line) at mid-rapidities (\y* \ < 0.35), using the px shift as Spx — {pt— < Pt >b 
) 15 /20, (solide line). The dashed line is obtained using Spx = (pt~ < Pt >t) 1 ' 5 /20 ioxpx < 7 GeV/c and 
Pt constant for px > 7 GeV/c. Right: Values of R^a u {pt) for minimum bias collsions at mid-rapidities 
{\y* \ < 0.35), using the p T shift as Sp T = (p T ~ < Pt > 6 ) 1 ' 5 / 2 0- The data are from PHENIX^l. 

We turn next to minimum-bias dAu collisions at central rapidity. Here < px >= 0.39 GeV/c. 
With the same value of n as above (n = 9.99), this corresponds to po = 1.346. Calculating the 
ratio dAu to pp we obtain the result of Fig. 1. At forward rapidity, RdAu decreases as y increases 
due two effects: 

The first effect is basically due to energy-momentum conservation. It has been known for a long 
time in hadron- nucleus collisions at SPS energies -low px "triangle"-. Its extreme form occurs in 
the hadron fragmentation region, where the yield of secondaries in collisions off a heavy nucleus 
is smaller than the corresponding yield in hadron-proton. This phenomenon is known as nuclear 
attenuation. It turns out, that, at RHIC energies, this effect produces a decrease of RdAu of 
about 30 % between rf = and rf = 3.2. 

The second effect is the increase of the shadowing corrections in dAu with increasing y* . This 
produces a decrease of RdAuipT) between rf = and rf = 3.2 of about 30 % for pions produced 
in minimum bias collisions. Therefore, we expect a suppression factor of about 1.7 between 
R(IAu{pt) at r]* = and at rj* = 3.2, practically independent of px- This is consistent with the 
BRAHMS results. 
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